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Abstract 

A form-finding problem for tensegrity structures is studied; given an abstract graph, 
we show an algorithm to provide a necessary condition for it to be the underlying graph 
of a tensegrity in M. d (typically d — 2, 3) with vertices in general position. Furthermore, 
for a certain class of graphs our algorithm allows to obtain necessary and sufficient 
conditions on the relative position of the vertices in order to underlie a tensegrity, for 
what we propose both a geometric and a symbolic approach. 

Subject classification: 05C85, Graph algorithms. 
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1 Introduction 

In this paper we study an instance of the so-called form-finding problems for tensegrity 
structures. These have brought special attention both among mathematicians and engineers 
since the seminal works of Kenneth Snelson around 1948 (see ^Hj). Roughly speaking, a 
form-finding problem for a tensegrity structure asks to determine a geometric configuration 
of points and straight edges in R d (typically d = 2,3) such that the whole structure is in 
a self-tensional equilibrium. The word tensegrity was coined from tension and integrity by 
Buckminster Fuller, deeply impressed by Snelson's work. 

Apart from a purely mathematical interest E3 , understanding these structures has 
applications to architecture and structural engineering JH] and has led to interesting models 
for viruses and cellular structures Cj • It is also considered a useful tool for the study 
of deployable structures jlUl 1141 117j . Previous works have proposed a number of different 
approaches to solve form-finding problems, which can be found in the recent review jllj . 

In particular, the present paper deals with the form-finding problem of building tenseg- 
rity structures with a given underlying graph G, in a given M. d . The graph has to be 
understood as an abstract graph, i.e., a set of vertices and pairs of vertices (edges). We aim 
to solve the following two problems: 

• First, to decide whether G can be the underlying graph of a tensegrity structure in W 1 . 



^In memoriam. 

'Departamento de Matematicas, Universidad de Alcala, e-mail: david.orden@uah.es Research partially 
supported by grants MEC MTM2005-08618-C02-02 and CAM S-0505/DPI/000235. 
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• In case such a tensegrity with underlying graph G is possible, to characterize the 
relative position of its vertices. 

In order to solve these problems, we first look for decompositions of tensegrities into basic 
instances, called atoms. This motivates a combinatorial method that allows to decompose 
a graph G into the smallest graphs that can underlie a tensegrity. In order to build up a 
tensegrity with graph G, we propose to reverse its decomposition: We show that this solves 
the above problems for a certain class of graphs and we present two different approaches. 
The first one looks at the geometric structure of the tensegrity; it is quite visual and provides 
intuition of the intrinsic properties of tensegrity structures. However, it becomes difficult to 
use for complicated structures. The second approach condenses in a matrix the information 
about the tensegrity; this allows to use tools from Symbolic Computation, despite being 
less intuitive. 

The paper is organized as follows: The basic notions and results are introduced in 
Section EJ Then, Section |3] introduces a method to decompose a tensegrity into atoms, 
which motivates a decomposition of the abstract graph G, reversed then by geometric 
means. Finally, in Section 0] a rigidity matrix is used for a symbolic resolution. 

2 Preliminaries 

In this section we introduce the basic notions and results used in the paper. Despite it 
aims to be self-contained, an interested reader can look at |18j for further examples and a 
more detailed overview of the mathematical concepts. Let us introduce first the rigorous 
definition of "self-tensional equilibrium" : 

Definition 2.1. Let G = (V,E) be an abstract graph: 

• A framework G(P) in M. d is an embedding of G on a finite point configuration P := 
{pi, . . . ,p n } m with straight edges. In the sequel we will focus on general position 
point configurations (no d+1 points lie on the same hyperplane). 

• A stress w on a framework is an assignment of scalars Wij (called tensions) to its 
edges. Observe that = Wji, since they refer to the same edge. 

• Such a w is called a self-stress if, in addition, the following equilibrium condition is 
fulfilled at every vertex: 



That is, for each vertex pi the scaled sum of incident vectors piPj is zero. 

Observe that the null stress is always a self-stress, of no interest for us. Note also that 
all scalar multiples of a self-stress (in particular its opposite) are self-stresses as well. We 
will see later that, indeed, the space of self-stresses on a given graph is a vector space. 




(1) 



ij edge 
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Lemma 2.2. Let p G P be a vertex of a d- dimensional framework G(P) such that P is in 
general position. Given a non-null self-stress on G{P), either at least d + 1 of the edges 
incident to p receive non-null tension, or all of them have null tension. 

Proof. The result is true for any d, but the reader may consider d = 2, 3 here. Let k be the 
number of edges incident to p that have non-null tension. The equilibrium condition on p 
implies having k vectors in R d , with common tail, which add up to the zero vector. For 
k < d+1, this is only possible if their k + 1 endpoints do not span a /c-space. But either 
k = or this contradicts the general position assumption. □ 

As a consequence, the next property makes particularly interesting the study of a certain 
family of general position frameworks, the so-called {d + l)-regular ones, for which a null 
tension on a single edge propagates to the rest of them: 

Corollary 2.3. Given a d- dimensional framework all of whose points are in general position 
and have exactly d+1 incident edges, a self-stress is non-null if, and only if, it is non-null 
on every edge. 

The following definition introduces our final object of study, which is a physical model 
of the mathematical objects defined above: 

Definition 2.4. We define a tensegrity structure T{P) to be a self-stressed framework in 
which: 

• Edges ij such that Wij > have been replaced by inextensible cables (its endpoints 
constrained not to get further apart), 

• Edges with < have been replaced by unshrinkable struts (endpoints constrained 
not to get closer together), and 

• Edges with = have been removed. 

If no confusion is possible, a tensegrity structure T(P) will be denoted by just T. For 
another physical interpretation, one can think of cables and struts as springs endowed with 
a certain tension, respectively inwards and outwards. That is; cables and struts incident to 
point pi have respectively tensions in the direction of +PiPj (inwards) and -piPj (outwards), 
see Figure [I] 




Figure 1: Left: two cables (+) and a strut (— ) incident to pi. Right: Their representation 
as springs with inwards and outwards tensions. 
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Observe that, given a tensegrity structure, it might be possible to replace the struts 
by bars which react to the surrounding tensions. For example, if we replace the strut in 
Figure H by a bar, this will receive an outwards tension at pi, as a reaction to the sum of 
cable tensions. Such a replacement is usual when constructing tensegrity sculptures, like 
those in [TH] . 

The most emblematic tensegrity structure is shown in Figure 121 Named oblique tri- 
angular prism with rotational symmetry, it is composed of nine cables, six of which form 
two copies of an equilateral triangle, the top one rotated 30 degrees, joined by three struts 
alternating the rest of cables. Thick edges denote the struts, which could be replaced by 
bars as before. 



The last definition in this section introduces the smallest tensegrities possible, which we 
will show in Section |3] to be the fundamental bricks for building up any tensegrity: 

Definition 2.5. We define a self-stressed atom in M. d (d = 2, 3) to be a general position 
realization of the complete graph JQ+2 , together with its unique (up to constant multiplica- 
tion) non-null self-stress. The tensegrity atoms are then obtained replacing edges by cables 
and struts. When no confusion is possible, we will just refer to atoms. 

Figure shows half of the possible tensegrity atoms in W d for d = 2,3, where thick edges 
denote struts. The other half is obtained by interchanging cables and struts or, equivalently, 
by considering the opposite tensions. It is not difficult to check, using Lemma 12. 2| that 
configurations with fewer points or edges do not admit self-stresses apart from the null 
one. The non-trivial fact that the above frameworks do admit a unique (up to constants) 
non-null self-stress appears in |12| . where existence is proved by the following result. 

Proposition 2.6. Let Yl?=i ^%Vi = 0; Y17=i Aj = £>e an affine dependence on a point set 
P = {pi, . . . ,p n }. Then, Wij : = AjAj defines a self-stress on the complete graph K(P). 



Pi 




Figure 2: The oblique triangular prism. 
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Figure 3: Types of tensegrity atoms in M? (left) and M 3 (right). 



Proof. For any ft G P, we have: 

n n n 

^2 Wij(pi - Pj) = ^ X i X j(Pi ~ Pj) = X iPi ^2 A J ~ Xi ^2 X i P i' 
ij£K j=l j=l j=l 

which equals zero. □ 

In order to prove the uniqueness up to constants, consider two different self-stresses, 
one not a scalar multiple of the other. Then some linear combination of them would cancel 
the tension at a particular edge but not at all of them, in contradiction with Corollary 12.31 
Note that we are using the claimed fact that self-stresses form a vector space, as will be 
shown at the beginning of the next subsection. 



3 Geometric approach 

In this section we present a geometric algorithm to decompose a tensegrity into atoms. This 
decomposition motivates a combinatorial one, which opens the way towards the resolution 
of the two problems posed in Section ^ 

3.1 Decomposing into atoms 

Let G(P) and G'(P') be two self-stressed frameworks such that PUP' is a point configuration 
in general position. Let w and w' be their self-stresses. For the framework G(P) U G'(P') 
obtained by union of vertices and edges, one can define the sum of self-stresses w + w' in the 
natural way: Assign tension Wij + w'^ to common edges ij and maintain the initial tension 
at the others. 

It is easy to observe that equations (^Q) are fulfilled and hence w + w' is indeed a self- 
stress. Furthermore, the space of self-stresses on a given graph G = (V, E) together with 
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this sum and the product by a scalar form a vector subspace of M E , when the latter is 
identified with the space of all self-stresses. 

Abusing notation, we denote by G + G' the self-stressed framework obtained. Observe 
that, after this addition is performed, one can appropriately replace edges by cables and 
struts in order to obtain a tensegrity structure T + T'. Hence, the sum of tensegrities yields 
another tensegrity. 

Observation 3.1. We will only consider this kind of addition when P and P' have at 
least d points in common; otherwise we obtain either two separate tensegrity structures or 
one of them hanging from the other. 

The main result in this section states that, reciprocally, under our conditions every 
tensegrity can be decomposed into a sum of tensegrity atoms: 

Theorem 3.2. (Atomic decomposition of tensegrities) Every non-null tensegrity 
structure T[P), P in general position, is a finite sum of tensegrity atoms. This decom- 
position is not unique in general. 

Proof. Let G{P) and w be the framework and non-null self-stress associated to T(P). We 
show how to obtain, by addition of atoms, a chain of non-null self-stresses w' on G(P) in 
which the number of vertices with only null incident tensions (null vertices) is increased at 
each step. At the end we come up with a self-stressed framework with only null vertices, 
so that the original tensegrity T will be the sum of the opposites of those atoms that have 
appeared in the process. 

Let us focus on the two-dimensional case, since the ci-dimensional one is carried out 
analogously: At each step, an arbitrary non-null vertex a £ P is chosen to be converted in 
a null one. By Lemma 12.21 only the following two cases are possible (see Figures 0] and |SJ) : 

• Type 1: If exactly three incident edges ab,ac,ad have non-null tension, we consider 
the atom K of vertices a, b, c, d. Since this atom has a non-null self-stress w K which is 
unique up to constants, we can choose w^ b to be the opposite of the tension assigned 
to edge ab at the current stress w', i.e. w^ b '■= —w' ab - Because of the equilibrium at a, 
it turns out that also w^ c = —w' ac and w^ d = —w' ad . Therefore, adding w to the 
current self-stress makes vertex a have only null tensions at incident edges (i.e. makes 
it disappear from the induced tensegrity) . See Figure EJ where dashed interior edges 
in the second picture are opposite to those in the first one. Note that at b, c, d the 
edges be, bd and cd may have appeared with non-null tension, but these extra edges 
do not affect a. However, we will be concerned about them later. 

• Type 2: If a £ P has incidence degree greater than 3, let b, c, d be neighbors of a. 
Consider the atom K of vertices a, b, c, d (and all the possible edges between them) 
and choose it to have tension w^ b := — w' ab at edge ab. Hence, obviously w' + w 
has null tension at edge ab. Again, other edges be, bd, cd may appear with non-null 
tension, but not incident to a. Hence, repeating this process if needed, we obtain a 
self-stress on G(P) in which a has only three incident edges with non-null tension (i.e. 
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Figure 4: Type 1 step, exactly three incident edges with non-null tension. 

in the induced tensegrity, a has only three incident edges). Now we are in the previous 
case. See Figure El where now only dashed edge ab is guaranteed to be opposite to its 
filled counterpart. 




Figure 5: Type 2 step, more than three incident edges with non-null tension. 

Since a sum of self-stresses is another self-stress, after a finite number of these steps we 
get a self-stress with at least one more null vertex, for which we can iterate the process 
until all vertices become null. Note that different choices of vertices to make null may lead 
to different decompositions. □ 

Remark 3.3. The reader should notice that the addition of an atom only changes the value 
of the self-stress on edges of the atom. Hence, adding an atom might cancel other tensions 
than the intended ones, but only at edges contained in the atom. 

Motivated by the geometric process in Theorem 13.21 we define now the following com- 
binatorial algorithm, that can be applied to any abstract graph G: 

Algorithm 3.4. (Combinatorial decomposition) 

INPUT: abstract graph G = (V,E) and dimension d. 

OUTPUT: list L of "atoms", where each atom is a subset of (d + 2) elements ofV. 

1. Initialize L = 0. 

2. While E is not empty, choose a vertex a £ V with minimum degree and: 
2.1 If a has degree < d, remove its incident edges from E. 
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2.2 If a has degree d+1, let ao, • • • , a^ be its neighbors. Remove the edges aai from E. 
Add to E all the edges aiOj that were not in E. Insert the atom {a, ao, . . . , a^} 
to the list L. 

2.3 If a has degree at least d + 2 do the following until it has degree d + 1, then go 
to 2.2: Choose d + 1 neighbors ao, . . . , a<2 of a. Remove the edge aao from E and 
insert to E all the edges aiaj that were not in E. Insert the atom {a, oq, . . . , a^} 
to the list L. 

3. Return L. 

See Figures |H1 and [7| for examples. This combinatorial algorithm is the tool for the first 
result towards the resolution of the problems posed in Section 1: 

Theorem 3.5. Given an abstract graph G, in order for it to underlie a tensegrity T(P) 
in M. d it is a necessary condition that, chosen a combinatorial decomposition, for every edge 
pq of G there is an atom containing both endpoints p and q. (In other words, that step 2.1 
of the algorithm removes only edges that were contained in the complete graph defined by 
some atom of the combinatorial decomposition). 

Proof. Choose a combinatorial decomposition as above. If such a tensegrity T(P) exists, 
then the combinatorial decomposition induces a geometric one: On the one hand, the 
geometric counterpart of step 2.1 shows up naturally because of Lemma 12.21 On the other 
hand, steps 2.2 and 2.3 correspond to the geometric steps of types 1 and 2, respectively, 
in the proof of Theorem 13.21 with tensions determined by the edge(s) to be deleted. The 
necessary condition in the statement is then a consequence of Remark 13.31 □ 

Figure El shows a combinatorial decomposition that fulfills the condition of Theorem 13 .51 
while Figure [7| shows one for which edge pq shows that the condition is not fulfilled, hence 
the graph cannot underly a tensegrity. 




Figure 6: Combinatorial decomposition fulfilling the condition of Theorem 13.51 

Therefore, when asked about the existence of a tensegrity in M. d with underlying graph 
G, we first get a combinatorial decomposition and then check the condition of Theorem 13 .51 
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Figure 7: Combinatorial decomposition not fulfilling the condition of Theorem 13.51 

If this is not fulfilled, then the answer is negative. However, in the next subsection we 
show that if the graph G admits a combinatorial decomposition of a certain type, then 
this decomposition can be geometrically reversed, leading in a straight way to a complete 
characterization of the point sets P on which G(P) admits a tensegrity. 

3.2 Geometrically solving form-finding problems 

Let us start defining a class of combinatorial decompositions, for which the reader can find 
an example in Figure 

Definition 3.6. We call edge-inserting combinatorial decompositions those in which all 
extracted atoms but the last one introduce at least one edge to the intermediate graph. 

The following result states that, for this class of combinatorial decompositions, their 
geometric reversal provides a solution to the problems posed in Section ^ Note the request 
of a self-stress non-null on every edge, in order for G to underlie the tensegrity: 

Theorem 3.7. If an abstract graph G admits an edge-inserting combinatorial decompo- 
sition, then the reconstruction of the graph from the atomic decomposition produces a set 
of equations and negated equations characterizing the point sets P in general position that 
make G(P) admit a self-stress non-null on every edge (and hence the ones that admit a 
tensegrity T(P) ). 

Proof. In order to determine which choices of coordinates underlie tensegrities, we have to 
consider as variables the coordinates of the new points pi added in the reconstruction (the 
first d + 2 points can be arbitrarily chosen, since tensegrities are projectively invariant |13j). 
Thus, the tensions of the edges in each atom introduced are functions on these variables. 

Furthermore, recall from Definition 12.51 that the self-stress of each atom has one degree 
of freedom. Hence, for the reconstruction of a general combinatorial decomposition we have 
to consider one extra variable ay for each atom (except for the first one, whose stress can 
be considered a normalization constant). 

However, these extra variables are not needed for edge- inserting combinatorial decom- 
positions: For each step of the geometric reconstruction, at least one edge was inserted by 
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the combinatorial decomposition algorithm and has to be removed in this precise recon- 
struction step. Then, the relative self-stress given to this atom (considered as a function 
on the positions) can be determined at the insertion step: It is exactly the one that cancels 
the tension(s) at the edge(s) to be removed. 

Therefore, the reconstruction in this case leads to a system of equations f(p\, ■ ■ ■ ,p n ) = 
and negated equations f(pi, ■ ■ ■ ,p n ) 7^ 0, one equation for each edge not in G that appears 
in an intermediate step of the process (expressing that this edge has tension zero and does 
not appear in the tensegrity) and one negated equation for each edge in G (expressing that 
it does appear in the tensegrity). □ 

Let us point up that two decompositions of G lead to equivalent collections of conditions, 
since the tensegrities constructed using one decomposition can always be decomposed and 
reconstructed using the other. Hence, they have to fulfill both sets of necessary and sufficient 
conditions, and therefore these have to be equivalent. In order to finish the section, we 
illustrate the decomposition and reconstruction process with the following example: 

Example 3.8. The graph G = ({1, 6}, {12, 14, 16, 23, 26, 34, 35, 45, 56}) underlies a 
tensegrity T(P) in general position in M 2 if, and only if, the triangles P1P2P& and pzP^P'o 
are in perspective position, i.e., the lines P2P3,P5P6 andp^pl are concurrent. 

We start with the edge-inserting combinatorial decomposition of G depicted in Figure El 
In Figure |H1 we show the reversal of this decomposition, where dashed edges are the inserted 
ones that have to be removed. 




Pi Pi Pi 



Figure 8: Reconstruction of an edge-inserting combinatorial decomposition. 



As observed in the proof of Theorem 13.71 the initial points P2,P4,,P5,P6 can be chosen 
arbitrarily. Then we add the atom P2,P3,P4,P5, whose self-stress is determined by the 
condition of canceling the tension at edge P2P5- The feasible positions of p% are characterized 
by the equations and negated equations corresponding respectively to the edges omitted and 
depicted in the middle picture. 

Finally, the addition of the atom pi,P2,P4,P6 has to remove edges and PaPq- This 
fact allows us to rephrase geometrically the set of equations and negated equations obtained 
for the graph G: The point p\ has to be placed on the line spanned by the resultant of the 
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"undesired" tensions W24 and uj^q so that the tension on edge pip^ replaces them. Moving 
pi along the feasible line towards the intersection of and pipi forces tension W26 at 
point P2 to be very small. But this equals tension wq2 at point p$; therefore, p\ has to be 
close to the intersection of pTPi and pEpe as well. 

Taking this argument to the limit, we conclude that p\ has to be chosen in such a way 
that the lines P2p3,pEPe and pTPi are concurrent. 



4 Symbolic approach 

This section is devoted to show that the form-finding problems considered can be reformu- 
lated in a matricial form, and how this can be used to obtain the equations and negated 
equations in Theorem 13 . 71 bv means of symbolic computations. In particular, this provides 
an alternative to the geometric ad hoc reasonings like the one in Example 13. 8( and allows 
to deal with more complicated and less visual examples. 



4.1 The rigidity matrix 

Before starting this subsection, let us point up that in this paper we are not dealing with 
the rigidity of tensegrity structures; we are just concerned by their self-equilibrium. An 
interested reader can see ^H] for more information on rigidity of frameworks. However, the 
following notion from rigidity analysis turns out to be useful for our purposes: 

Definition 4.1. Let G(P) be a framework with n vertices and e edges in M. d . Its rigidity 
matrix R(P) has e rows and nd columns, defined as follows: 

• There is a row per edge ij of the framework, with % < j and in lexicographic order. 

• Each block of d columns is associated to a vertex pi and contains either the d coordi- 
nates pi — Pj, at those rows corresponding to edges ij incident to pi, or zeros at the 
rest of rows. 

For a complete framework on n vertices (with all the possible edges), the rigidity matrix 
has the following condensed form: 



vertex pi 

I ' 

( P1-P2 Pi-pi 

Pi ~P3 

V 





P3-P1 



edge 12 
edge 13 



Pn-l-Pn Pn-Pr, 



J 



The key observation is that the equilibrium equations Q can be restated in matricial 
form as 

w ■ R(P) = (2) 

where w is a 1 x e vector of entries Wij,i < j (recall that Wij = Wji) and the right-hand 
side is the 1 x nd zero vector. That is to say, self-stresses w are row dependencies for the 
rigidity matrix, what leads to the following observation: 
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Observation 4.2. For a framework G(P), being the underlying graph of a tensegrity T(P) 
is equivalent to the existence of a w with no null component in the left kernel of the rigidity 
matrix R(P) of the framework. 

This observation already suggests a method to characterize the point sets P in general 
position that make G(P) admit a tensegrity T{P): Computing the kernel of the rigidity 
matrix R(P) for the framework G{P). For the sake of consistency, we refer the interested 
reader to for an example of the use of this method, and we show here instead how 
to use the decomposition-reconstruction method together with the matricial expression of 
equilibrium. 

4.2 Symbolically solving form-finding problems 

According to Observation 13. 11 the reconstruction of a combinatorial decomposition inserts 
at most two new points at each step, with indeterminate coordinates x, y 6 M. d . For the 
framework to maintain equilibrium after their addition, the points have to belong to the 
projection of the solutions of (J2j) over the space of the variables Xi,yj. Furthermore, in order 
to look for non-null tensions on every edge, the inequations uiij ^ for every i,j have to 
be considered. Equivalently, extra variables tij and equations lo^ty = 1 can be introduced 
to the system. 

It is well-known (see Lemma 1 in page 120 of jl]) that the polynomial elimination of 
variables Wij,tij from the above system of equations contains the above projection. There- 
fore, polynomial elimination provides necessary conditions for the position of the vertices, 
which can be later symbolically tested for sufficiency. 

We illustrate this symbolic method with the following example, in which in order to 
perform polynomial elimination we have chosen to compute the intersection of M\x, y] with a 
Groebner basis of the equations, for an elimination order in which variables Xi,yj are smaller 
than the Wij,tij's; see Theorem 2 in page 113 of jl]. Although there are more efficient 
elimination techniques, the reason for this choice is to present the operations performed 
with the software Maple \9'. which, apart from being broadly known, has a specific linalg 
package for linear algebra. For more efficient elimination software, the reader can consider 
CoCoA 2 . Let us also point up that the graph considered has eight triangles and, in fact, 
is the same as the graph of an octahedron. 

Example 4.3. The graph G = ({1, . . . , 6}, {12, 13, 14, 15, 23, 25, 26, 34, 36, 45, 46, 56}) of an 
oblique triangular prism underlies a tensegrity if, and only if, the six vertices lie on a ruled 
hyperboloid that contains the edges of one of the three cycles of length 4 in G. Equivalently, 
if and only if the planes containing four alternating triangles intersect. 

We start with the combinatorial decomposition of G in which the vertex 6 is chosen at 
step 2 of Algorithm 13.41 This gives the atoms 2,3,4,5,6 and 1,2,3,4,5 (see Figure EJ). In 
the sequel we detail the steps of the reconstruction process using Maple 7, and omitting 
some outputs of no interest: 

• According to the proof of Theorem 13. 71 the first d + 2 points pi,...,P5 can be arbitrarily 
chosen: 
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> pl:=[0,0,0] : P 2:=[l,l,l]: p3:=[0,l,0]: p4:=[l,0,0]: p5:=[0,0,l] 



• Solving the corresponding equation (j2J we get the tensions of the atom pi,P2,P3,Pi,P5- 
In order to generate the equations, we use the command geneqns of the linalg package, 
which needs the equivalent transpose form jR* • w t = of iffi) . Then we solve them and, 
since according to the proof of Theorem 13.71 the stress of the first atom can be considered 
a normalization constant, we take the value 1 for the parameter obtained. 

> d:=3: zeros : = [0 , , 0] : with (linalg) : 

> R:=matrix(10,5*d, [ 

op(pl-p2) ,op(p2-pl) ,op(zeros) ,op(zeros) ,op(zeros) , 
op(pl-p3) ,op(zeros) ,op(p3-pl) ,op(zeros) ,op(zeros) , 
op(pl-p4) ,op(zeros) ,op(zeros) ,op(p4-pl) ,op(zeros) , 
op(pl-p5) ,op(zeros) ,op(zeros) ,op(zeros) ,op(p5-pl) , 
op(zeros) ,op(p2-p3) ,op(p3-p2) ,op(zeros) ,op(zeros) , 
op(zeros) ,op(p2-p4) ,op(zeros) ,op(p4-p2) ,op(zeros) , 
op(zeros) ,op(p2-p5) ,op(zeros) ,op(zeros) ,op(p5-p2) , 
op(zeros) ,op(zeros) ,op(p3-p4) ,op(p4-p3) ,op(zeros) , 
op(zeros) ,op(zeros) ,op(p3-p5) ,op(zeros) ,op(p5-p3) , 
op(zeros) ,op(zeros) ,op(zeros) ,op(p4-p5) ,op(p5-p4) ] ) : 

> Rt :=transpose(R) : 

> eqs : =geneqns (Rt , [wl2,wl3,wl4,wl5, w23 , w25 , w26 , w34 , w36 , w45 , w46 , w56] , 
vector (n*d,0) ) : 

> tensions : =solve (eqs , wl2 , wl3 , wl4 , wl5 , w23 , w24 , w25 , w34 , w35 , w45) ; 

tensions := {wl5 = — 2u45,u;35 = wA5,wl2 = 2w45,wlA = —2w45,w23 = — w;45, 
u;34 = w45, u;13 = -2 w45, w45 = w45, w25 = -w45, w24 = -w45} 

> tensions : =subs (w45=l , tensions) ; 

tensions := {1 = l,tu35 = l,w\2 = 2,wU = -2,w23 = -l,wU = l,u>13 = -2, 
w25 = -l,u;24 = -1,^15 = -2} 

• Then we have to add the atomp2,P3,£>4jP5>P6) i n which pq has unknown coordinates x,y,z 
and the stress of the atom is determined by the cancelation of tensions at edges P2P4 and 
P3pE- The same operations as above lead to a second system of equations eqs2, in which we 
substitute -u;24 = l,u>35 = — 1 to get 

eqs2subs := {w22, + u>25 — w2Q x + w2Q = 0, — w23 — u>34 — x w36 = 0, w3A + ^45 — w4G x + 
w46 = 0, — w25 — wA5 — x w56 = 0, — w26 + w26 x + x w36 — wAQ + wA&x + xwb& = 0, — w26 + 
w2Q y — u>36 + t«36 y + y w46 + y w56 = 0, — u>26 + w2Q z + z w36 + z w46 — w56 + w56 z = 
0,1 + u>25 - w2§y + w2Q = 0,u;23 + 1 - w2§z + w26 = 0,ti;34 - 1 - w36y + w36 = 
0, -w23 + 1 - z w36 = 0, -1 - it;34 - y u>46 = 0, -1 - w45 - z w46 = 0, -w25 + l-y w56 = 
0, -1 + wA5 - w56 z + w56 = 0} 

• In order to look for non-null tension on every edge, Corollary 12 . 31 turns out to be crucial; if 
there is a non-null self-stress, then all its tensions are null. Therefore, the self-stress obtained 
when adding the two atoms is non-null since wl2 = 2 ^ and this tension is not affected 



13 



by the second atom. In particular, it is not needed to consider extra variables tif In order 
to obtain the polynomial elimination, we compute a Groebner basis of the polynomials in 
eqs2subs, which we call polys2subs, for an elimination order in which variables x,y,z are 
smaller than the w^'s. The command gbasis of the Groebner package is used: 

with (Groebner) : 

G:=gbasis(polys2subs,lexdeg( [w23,w25,w26,w34,w36,w45,w46,w56] , [x,y,z] )) ; 

G := [x 2 — x — z 2 — y 2 + y + z, . . . and 40 polynomials more, involving Wij's. 



We conclude that the hyperboloid x 2 



y 



x + y + z = contains the set of 



points pq := (x,y,z) such that the framework admits a self stress non-null on every edge. 
Therefore this is a necessary condition. In order to test its sufficiency, we have considered 
an extra edge in the framework and forced its tension to be null, obtaining the following 
stress from the left-kernel of the corresponding rigidity matrix: 



w 



I y + z — 1 — x 
y + x - z + 1 
y + x - z + 1 
y + x - z + 1 
— y + x 
x — z 

1 

y(-y+x-z+l) 

y+z+x—1 
_ —y+z+x—1 
y+z+x—1 
z(—y+x—z+l) 

y+z+x—1 
—y+x—z+1 
y+z+x—1 
y—z—l+x 
y+z+x—1 



V 



We omit the (easy) computations checking that: 



• Under the condition x 2 
self-stress. 



x + y + z = 0of the hyperboloid, this is indeed a 



• The denominator is not null: We observe that points ^3,^*4 and p§ already lie on the 
plane y + z + x — 1 = 0. Therefore, if point p§ lied on that plane, the configuration 
would not be in general position. 

• Similar arguments prove that all the components of w are non-null. 

In conclusion, points p$ = (x, y, z) for which a tensegrity exists are precisely those lying 
on x 2 — y 2 — z 2 — x + y + z = 0. We now observe that this hyperboloid contains the edges 
of the quadrilateral P2P3P4P5- Indeed, it is the only one passing through the initial five 
points and containing those four edges, since the space of all hyperboloids has dimension 9. 
In order to conclude the first condition in the statement we just have to observe that the 
combinatorial decompositions that at step 2 of Algorithm 13.41 choose vertices 1 or 6, lead 



14 



to the 4-cycle 2345, those choosing 3 or 5 lead to 1264, and those choosing 2 or 4 to 1365 
(see Figure HJ). 

For the equivalent condition in the statement, which appears in 0, easy computations 
check that the four planes defined by P1P2P3, P1P4P5, P3P4P6 and P2P5P6 intersect precisely 
for x 2 — y 2 — z 2 — x + y + z = 0, and the same condition is obtained for the other four 
alternate planes. Let us finally note that an interested reader can find in [H] a different 
symbolic approach to the resolution of this problem, that uses comprehensive Groebner 
basis. 

Remark 4.4. We have to call the attention of the reader to the fact that, in spite of dealing 
with more complicated problems than the geometric approach, the usefulness of the algebraic 
one is also limited by the size and type of the problem. For instance, if Corollarv \2.^ cannot 
be used and all the inequations Wij 7^ have to be considered in order to look for non-null 
tension on every edge. This would introduce more auxiliary variables £12,^13, • • • ,^56 and a 
polynomial (1^12*12 — 1)(^13*13 — 1) • • • (^56^56 — 1); making the computations infeasible. 

Finally, let us recall that the above symbolic computations follow the decomposition- 
reconstruction method from Section We already noted that following Observation \4- 6 A 
it is also possible to compute symbolically the kernel of the 12 x 18 rigidity matrix R{P) 
for the whole framework G{P) . Instead, with the decomposition-reconstruction approach we 
compute the kernel of the two 10 x 15 submatrices corresponding to the two atoms. The main 
advantage of the decomposition reconstruction method is related to this fact: Since w24 and 
w35 were no longer variables and the system had fewer equations, the computation of the 
Groebner basis was easier. In addition, we did not need to introduce an extra variable tij. 

Acknowledgements 

Both authors would like to thank the initial help and later encouraging comments and 
revisions by professors Juan Llovet and Tomas Recio, specially about the algebraic part. The 
second author also wants to thank an anonymous reviewer for several valuable comments 
that helped improving the results in the paper. 

Final note 

The present paper is dedicated to the memory of Miguel de Guzman, who left us before 
it was finished. This work started when he gave a conference about tensegrities at the 
Departamento de Matematicas of the Universidad de Alcala in November 2003. Five months 
later he died after a seminal version had been finished. Although part of the results presented 
here still had to be formulated and proved in their final form, their essence was already 
contained in that preliminary version. 
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